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Abstract. Application of the theoretical results of the author are presented for the 
case of meteor streams of comet Encke. Theoretical determination of meteor orbit for 
comet Encke, as a parent body, is presented. Four significant theoretical meteor streams, 
corresponding to Tauds N, Tauds S, /3 Tauds and ^ Perds are found. The meteor stream 
membership criterion is applied to the photographic orbits of The lAU Meteor Data 
Center in Lund: Taurid meteor stream is found for several possible areas in phase-space 
of orbital elements. 
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1. Introduction 

Application of the meteor stream membership criterion and the method of prediction of 
theoretical radiants, suggested in Klacka (1995a: Meteor stream membership criteria - 
paper I) and Klacka (1995b: Meteor streams and parent bodies - paper II), is presented 
for the case of meteor streams connected with the comet Encke. 

Gravitational forces acting on meteoroids of the comet Encke yield that semimajor 
axis a (energy E, E = 1/a - the sign is not important, now), z— component of the 
angular momentum ( L^ = H^ = ■\/a(l — e-^) cos i ) and the angle tt {tt — u; + fl) 
show small dispersions during evolution for several 10,000 years. This suggest that these 
quantities should be used in determining the meteor stream membership and in prediction 



of theoretical radiants, e - eccentricity was taken as the fourth important quantity for 
determination of meteor stream. The reason for this choice is that distributions in E, L^, 
■K and e may be approximated with normal distribution. The data for the photographic 
meteor orbits of The lAU Meteor Data Center in Lund were used. 

2. Theoretical prediction of radiants 

Let us consider that comet Encke is the only parent body for a meteor stream. We use 
the result of section 6 of paper I in the form of "Energy —Hz — tt method" (see section 5.1 
in paper II) for the purpose of finding of meteor stream(s). The local minima yield four 
theoretical radiants, summarized in Table 1; the theoretical radiants should correspond 
to (3 Tauds, Tauds S, Tauds N and ^ Perds. 

Table 1. Theoretical radiants for comet Encke as a parent body. Energy-H^-Tr method. 

q[Au] e z[°] Lon »["] an 5n x^n lAjfp 

0.326 0.853 6.6 242 278 

0.324 0.854 4.8 118 42 

0.321 0.855 298 222 

0.321 0.855 62 98 
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3. Photographic data and membership to meteor stream 

Section 5 of paper I is used in determining the membership to meteor stream. The data 
for the photographic meteor orbits of The lAU Meteor Data Center in Lund were used. 
We have found that orbital elements E^ Lz, tt and e may be approximated with normal 
distribution. Thus, the density function in Eq. (8) of paper I is taken in the form of four- 
dimensional normal density function. The area il in Eq. (8) of paper I may be chosen in 
an infinity ways - we take, in this paper, only some of special cases summarized by the 
form 

j^p^Y^(\Q^_Q2iAX <Di, (1) 

where Qi = E, Q2 = e, Qt, = Hz, Qa = tt, and, the value Dc corresponds to the 
case when multidimensional normal distribution yields that meteors with D > Dc are 
members of a group, for given r^ (i = 1 to 4) and p, only with probability less than 5 
% {QiM corresponds to the mean value of the quantity Qi). We consider the following 



cases, as examples: 

1. p = 2, ri = r2 = rs = ri = 1. 

2. p = 1, ri = r2 = rs = 1, r^ -^ oo, 

3. p = 2, ri = r2 = rs = 1, r^ -^ oo, 

(We must stress that Eq. (1) is principally different from Eq. (1) in Klacka (1992) - the 
case 1 corresponds, roughly, to Eq. (1) in Klacka (1992); however, even in the case that 
area fl is of the type of a quadratic form, it may be even of indefinite type.) 

3.1. Description of the procedure and basic results 

As we have already mentioned, the data for the photographic meteor orbits of The lAU 
Meteor Data Center in Lund were used. The requirement 0.15 AU < q < 0.5 AU and 0.4 
< e < 0.999 and i < 20° yields 418 objects of the total set. The set of 418 meteors was 
used in finding the required meteor stream. 

The case 1, already defined, was ellaborated in the following way. At first, 418/36 = 
11.6 objects per 10°. The interval of tt, in which number of objects per 10° is greater 
than 11.6, corresponds to tt e (110°, 190°) = / (no object exists for tt e (105°, 110°)). 
The distribution for tt G / may be approximated by normal. Thus, initial set of 418 
objects reduced to 213 objects with tt G /. Four-dimensional normal distribution was 
used (quantities Qi — E, Q2 — e, Q3 — Hzi Qi ~ t^)- Since only several objects had 
8.5° < i < 20° (no object existed for i G ( 8.5°, 9.5°)), we reduced the initial set of 
213 objects into a set with i < 8.5° - the final starting set contained 191 objects. Four- 
dimensional normal distribution was applied to the set of 191 objects (quantities Qi — E, 
Q2 = e, Q^ = H^, Q4 = tt). The final result is characterized by 4-dimensional normal 
distribution with mean values Qimi Q2Mi Q3M, Qam ^nd covariant matrix of the vector 
(Qi, Q2, Qzi Qi) with components (Jij {i, j = 1 to 4) ~ see Table 2. We have obtained 
130 objects for the case 1 (probability 95 % - Dc — 3.31). (Quantities ai, i — 1 to 4, used 
in Eq. (1) are given by the relations ai = y/oii, i = I to 4.) 

In the cases 2 and 3, already defined, the three-dimensional normal distribution was 
used (quantities Qi = E, Q2 ~ e, Qj, ~ H^). In both cases the statistics yielded that the 
obtained distributions in n corresponded to normal for n G (110°, 190°) = /, the number 
of objects with tt G /' was smaller than the number of objects with tt G /. As for the 
distribution in inclinations, the same holdcd as in the case 1. Thus, the final starting set 
of the data contained 191 objects. Three-dimensional normal distributions were applied 
to the set of 191 objects (quantities Qi = E, Q2 = e, Q^ = H^). The final results are 
characterized by 3-dimcnsional normal distributions with mean values Qim, Q2M, QaM 



Table 2. Mean values and components of covariant matrix for the case 1. (Probability 
95 %^Dc^ 3.31). 
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and covariant matrix of the vector (Qi, Q2, Q3) with components ct^ (i, j = I to 3) — 
see Table 3 for the case 2 and Table 4 for the case 3. We have obtained 130 objects for 
the case 2 (probability 95 % - Dc — 4.69) and 131 objects for the case 3 (probability 95 
% - Dc ^ 3.02). (Quantities ai, i — I to 3, used in Eq. (1) arc given by the relations 
(Ti = y/o^, i = 1 to 3.) 



Table 3. Mean values and components of covariant matrix for the case 2. (Probability 
95 %- Dc = 4.69). 



QlM 


Q2M 


Q3M 


0.4672623 


0.8323077 


0.8076754 


0-11 


0'12 


CTlS 


0.0034363 


-0.0006934 


-0.0014475 


CT21 


(722 


(723 


-0.0006934 


0.0007924 


-0.0011204 


(731 


(732 


(733 


-0.0014475 


-0.0011204 


0.0037191 



Table 4. Mean values and components of covariant matrix for the case 3. (Probability 
95 %^Dc^ 3.02). 
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3.2. Comparison of the results 

Table 5 presents objects classified as members of all three meteor streams - each object 
belongs to the stream in the cases 1, 2 and 3, simultaneously (127 objects). 

Table 6 yields those objects which are members only in one or two streams; the term 
"yes" denotes that an object is member of a given stream, the term "no" denotes that 
an object is not member of a given stream. 

While in the cases 1 and 2 the number of objects belonging only to one of the streams 
to the number of objects belonging to the both streams is 3/127 « 2.4%, all three cases 
yield similar ratio 7/127 ss 5.5%. It is evident that this number may even increase when 
other areas Q. are used. 

All three cases yield Aq e (165°, 265°). Distribution in Xq is not approximable by 
normal distribution. The relation between quantities of different distributions in a form 
of linear regression is not wise. This holds, as a special case, between an orbital element 
and Aq - even in a small interval AAq (as for larger interval, the procedure may even 
not be used due to the fact that periodicity is not secured). Evident approximate relation 
exists between tt and Aq due to the well-known relation among tt — il — Aq; however, 
we again emphasize that density functions for tt and Aq are different. However, in order 
to be able to compare our Taurid stream with that obtained in the past, we present the 
following linear regressions: 

q[AU] = (0.096 ±0.018) - (0.791 ± 0.811)10"''Aq[°] , (2) 

7r[°] = (-12.065 ±4.710) + (0.7563 ± O.O213)A0[°] , 



Table 5. Objects belonging to all three meteor streams (all three cases). 

2192 2198 2209 2254 2256 2259 2264 2372 2381 2391 

2418 2452 2457 2461 2480 2485 2506 2516 2526 2528 

2532 2552 2556 2565 2576 2580 2586 2604 2611 2629 

2637 2668 2680 2681 2684 2702 2713 2720 2722 2724 

2726 2733 2739 2747 2748 2752 2755 2761 2763 2764 

2768 2769 2777 2787 2793 2794 2796 2801 2805 2807 

2812 2814 2815 2816 2818 2820 2821 2822 2823 2826 

2834 2840 2842 2844 2846 2847 2849 2850 2852 2854 

2857 2858 2864 2866 2869 2874 2875 2877 2878 2886 

2890 2891 2892 2899 2900 2903 2906 2907 2910 2911 

2923 2935 2979 2980 2994 3000 3001 3012 3029 3036 

3041 3043 3050 3063 3066 3069 3082 3084 3093 3112 

3116 3125 3134 3136 3178 3189 3264 

Table 6. Objects of meteor streams maximally in two of the three cases. 
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for the case 1, 

q[AU] ^ (0.094±0.018) - (0.671 ± O.816)1O"''A0[°] , (3) 

7r[°] = (-18.410 ±5.246) + (0.7826 ± O.O237)A0[°] , 

for the case 2, and, 

q[AU] ^ (0.087±0.017) - (0.363 ± O.791)1O"^A0[°] , (4) 

7r[°] = (-16.756 ±4.887) + (0.7759 ± O.O221)A0[°] , 



for the case 3. Comparison of Eqs. (2) - (4) with the results of Porubcan and Stohl (1987 
- Table 1) shows that our Taurid meteor complex is different from that obtained on the 
basis of D-criterion of Southworth and Hawkins (1963). 

We must stress again that criterion suggested in paper I (in section 5), and applied 
in this paper, is of purely mathematical character, based on mathematical statistics. Our 
criterion is not of physical character, as was interpreted in Neslusan (1996). 

4. Conclusion 

We have applied methods developed in papers I and II to meteor streams of comet Encke. 
We have found, in a simple physical way, theoretical radiants which should correspond 
to comet Encke as a parent body. The coincidence with the reality is not very bad. 

Taurid meteor complex, based on the data of photographic orbits, was found for 
various areas in phase space. Comparison of the results based on statistical method of 
multidimensional density function is not consistent with that obtained on the basis of 
Southworth and Hawkins D-criterion. The difference will even increase with the increasing 
number of photographic meteors. 

We can conclude that our definition of the Taurid meteor complex corresponds to the 
multidimensional normal distribution in phase space of orbital elements. As an example, 
we have used three different areas fl. Physical problem is to make a decision which kind 
of area fi is the most reasonable. 
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